WEAK CONVERGENCE OF FINITE ELEMENT 
APPROXIMATIONS OF LINEAR STOCHASTIC EVOLUTION 
■^ ' EQUATIONS WITH ADDITIVE NOISE II. FULLY DISCRETE 

SCHEMES. 

o 

, MIHALY KOVACS, STIG LARSSON 1 - 2 , AND FREDRIK LINDGREN 1 

<D ' 

**■* • Abstract. We present an abstract framework for analyzing the weak error of 

g\ ' fully discrete approximation schemes for linear evolution equations driven by 

additive Gaussian noise. First, an abstract representation formula is derived 
for sufficiently smooth test functions. The formula is then applied to the wave 
equation, where the spatial approximation is done via the standard continu- 
^^ ous finite element method and the time discretization via an /-stable rational 

>— _ _ approximation to the exponential function. It is found that the rate of weak 

^— I i convergence is twice that of strong convergence. Furthermore, in contrast to 

pH ■ the parabolic case, higher order schemes in time, such as the Crank-Nicolson 

,j__j ' scheme, are worthwhile to use if the solution is not very regular. Finally we ap- 

ply the theory to parabolic equations and detail a weak error estimate for the 
linearized Cahn-Hilliard-Cook equation as well as comment on the stochastic 
heat equation. 

> 

£^ . 1. Introduction 

CN ■ Let U, H be real separable Hilbert spaces and consider the following abstract 

£f) ' stochastic Cauchy problem 

O ' 

£s| . (1.1) dX(t)+AX{t)dt = BdW(t), i > 0; X(0)=X , 

where —A is the generator of a strongly continuous semigroup {E(t)}t>o on ^> 
B € B(U, %), where B(U, TL) denotes the space of bounded linear operators from U 
to T-L. The process {W(t)} t >o is a U- valued Wiener process with covariance operator 
5_i , Q > (selfadjoint, positive semidefinite) with respect to a filtration {Tt}t>o on a 

probability space (£2, T, P). We note that, strictly speaking, the process {W(t)}t>a 
is ^/-valued if and only if TrQ < oo. Finally, Xq is an J^-measurable H-valued 
random variable with finite mean. If 

(1.2) Tr (/ E(t)BQB*E(t)*dt\ < oo, 
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then the unique weak solution is given by (see [HJ Chapter 5]) 

(1.3) X(t) = E(t)X + I E(t- s)B dW{s). 

Jo 

With G: "H — > K being a twice Frechet difFerentiable function with bounded and 
continuous first and second derivatives, we study the weak error 

(1.4) e(T)=E(G(X(T))-G(X(T))), 

where X(T) is some approximation of the process X at time T. 

This paper is a sequel to [15]. It consists of three parts. In Section [2] we define 
some central concepts and state important background results used throughout the 
paper. In Section [3] we show that the error formula in [151 Theorem 3.1] holds 
for a much wider class of approximations of the solution to (|1.1[) than stated in 
that paper, which is concerned with spatial semidiscretization by finite elements. 
Finally, in Sections |4] and El the usefulness of the general error formula in Theorem 
13-ll is demonstrated through the fact that it can be applied to analyze fully discrete 
schemes for a wide class of stochastic evolution equations: hyperbolic and parabolic 
alike. The basic line of proof of the general formula is adapted from [10] which is 
concerned with the stochastic heat equation. 

The statement of Theorem [3J] deserves a motivation. Consider, for example, the 
case when X(T) is the value of a semidiscretization in space with finite elements, 
so that 

X(T) = E{T)X + I E(T - s)B dW(s) 
Jo 
has the same form as the solution (jl.31) of the original problem (jl.lj) . Then an error 
formula may be derived with the aid of Kolmogorov's backward equation and Ito's 
formula as in [15] . It turns out that the error analysis is substantially simplified if 
new processes are constructed by multiplying X{t) and X{t) by suitable integrating 
factors. That is, define new, drift-free processes 

E{T)X + I E(T-s)BdW{s) 
Jo 



(1.5) 


Y(t)- 


= E(T - t)X{t) 


and 






(1.6) 


Y(t)-- 


-- E(T - t)X{t) 



E(T)X + J E(T-s)BdW{s) 
Jo 



(1.7) X(T) = Y(T), X(T) = Y(T), 
with {Y(t)} t >o being the solution of the equation 

(1.8) dY(t)=E(T-t)BdW(t),t>0; Y(0) = E{T)X 
and {Y(t)}t>o solves 

(1.9) dY(t) = E(T - t)B dW(t), t > 0; f(0) = E(T)X Q . 

However, if X{T) is the result of a time-stepping scheme, then a process of the form 
(jl.3[) is not immediate. We note that in [S] the interpolation of the time-stepping 
operator family is performed in a manner that results in a family of deterministic 
operators {E(t)}t>o with a weaker type of semigroup property. This property is 
sufficient to mimic the computations in (|1.6|) . but the operator family does not 
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naturally admit deterministic error estimates. Thus, we have chosen an alternative 
path, as in [10J , where the discrete semigroup property of the time-stepping operator 
family is used first, followed by piecewise constant interpolation between the grid 
points. This will, as we shall see, yield a drift-free process as in the right hand 
side of (ll.6|) such that (|1.7[) still holds and known deterministic error estimates 
can be used almost immediately. In this case {E(t)} t >o will not be continuous in 
t and does not have the semigroup property but it turns out that these are not 
necessary. All we need to obtain the fundamental formulas (|3.6|) — (|3.8|) for the error 
is to assume that there exists a well-defined process {^(£)}*e[o,T] of the form 

(1.10) Y(t) = E(T)X + J E(T-s)BdW(s) 

Jo 

such that (TOl) holds. Here {E(t)} te[ o, T } C B{S,S) and B £ B(U,S), where S is 
a Hilbert subspace of H with the same norm (typically S = T~L or S is a finite- 
dimensional subspace of T-L). The process Y is then well defined if 

(1.11) Trf/ E(t)BQ[E(t)B}* dt) < oo, 

and in this case it will be the unique weak solution of (|1.9[) . Here the adjoint 
[E(t)B]* = B*E(t)* : S — > U is taken with respect to the scalar product in H. 

In Section 21 we apply the general formula from Theorem l3.1l to the error analysis 
of semi- and fully discrete numerical schemes for the stochastic wave equation 

(1.12) dU(t) - AU(t) dt = dW(t), U{t)\ dv = 0, t > 0; U(0) = U , J7(0) = V , 

where the solution process {U(t)}t>o and the Wiener process {W(t)} t >o take values 
in U = L2(T>), where T> is a sufficiently nice open bounded domain in W. d . Writing 
X(t) = [X 1 (t),X 2 (t)f := [U(t) : U(t)] T , X = [X . l7 X . 2 ] T ■- [U (h V ] T and 
A := —A, the wave equation (|1.12[) can be written in the form (|l.ip with 



A:= 



-I 

A 



B := 



It is well known that —A is the generator of a unitary, strongly continuous semigroup 
(and thus a group) on the space H = L 2 (T>) x (Hq(D))* . 

The first result in Subsection [4J] is a bound of the error e(T) in the more specific 
context of the wave equation but keeping the approximating process general. The 
bound is expressed in terms of the operators and initial data in (jl.lOp and (|1 .3[) . In 
Subsection 14.21 we apply this to single step rational approximations of (jl.ip , that 
is, to solutions of the scheme 

X ] = R(kA)(X : >- 1 +B(W(t ] ) - Wfa-i))), 

where k is the step size and where the rational function R fulfills the approximation 
and stability properties 

my) - e-*| < C\y\v+\ \y\<b, 
\R(iy)\<l, yGR, 

for some positive integer p and some b > 0. If the initial value is smooth enough 
and 

(1.14) HA^-^QA-^HTr < oo, 
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where || • ||tt denotes the trace norm, then for the weak error in (|1.4|) we have 

\e(T)\ = 0(k min( ^ 2m ) asfc^O. 

It is important to note that, in contrast to the stochastic heat equation, the conver- 
gence rate is improved with higher order schemes if the noise is sufficiently irregular 
because -^j increases with the order, p, of the method. This feature appears for 



fully discrete schemes investigated in Subsection 14.31 as well, where the spatial ap- 
proximation is performed with standard, continuous, piecewise polynomial finite 
elements of order r and the temporal approximation again with rational functions 
as in (|1.13p . If h denotes the size of the finite element mesh and if f|l .14[) holds, 
then for the first component X\ — U we have 

|E(G(Xi(T)) - G{X X {T)))\ = 0(fc m!n( ^ 2/M) ) + 0(h min{ ^+i 2M ) as M ^ 0. 

It is a general phenomenon that for non-smooth noise, the rate of weak conver- 
gence is twice that of the strong (mean-square) convergence. To be able to compare 
the weak rate with the strong rate of fully discrete schemes for the wave equation, 
and also because strong results for such schemes are absent in the literature, we 
prove a strong convergence result in Section [4.41 for the same algorithm. Indeed, we 
find that the strong rate is half the weak rate also in the case of the wave equation. 

Our main motivation in applying the results of Section [3] to the stochastic wave 
equation is the fact that this is one of the canonical SPDE's that is less understood 
from the numerical point of view. Results on weak convergence can be found in 
[13] but are confined to a leap-frog scheme in one spatial variable. Also, the test- 
functions in that paper differ from ours, making the results difficult to compare. In 
[15] the stochastic wave equation is discretized with finite elements in space only 
but in several spatial dimensions. The findings are in accordance with the results 
of the present paper. 

For the stochastic wave equation on the one-dimensional real line with white 
noise the strong rate 1/2 for the leap-frog scheme is proved in [25]. This rate is also 
proved to be optimal. For white noise in one dimension the algorithm discussed in 
the present paper only gives a strong rate of p/2(p+ 1) and hence it is not optimal. 
The reason for this is explained in [16] with the fact that the Green's function of 
the leap-frog scheme coincides at the meshpoints with the Green's function for the 
wave equation, which is not true for the present scheme. The latter paper, [16] . 
studies spatial discretization with finite elements in several dimensions with findings 
in agreement (in one dimension) with the finite difference spatial approximation 
studied in 23^ 

In connection to hyperbolic equations [5] should be mentioned, where the authors 
are concerned with weak as well as strong convergence of a time discretization 
scheme for a nonlinear stochastic Schrodinger equation. 

Finally, in Section [5j we apply Theorem 13.11 to the backward Euler in time and 
finite element in space approximation of the linearized Cahn-Hilliard-Cook equa- 
tion. The reason for doing this is twofold. First, we demonstrate that the general 
error representation formula is useful in the parabolic setting as well. Second, the 
stochastic heat equation, which is the canonical parabolic equation, has been stud- 
ied in 10 . While our general approach would certainly be applicable to the heat 
equation as well, it would just reprove a known result, maybe with a more trans- 
parent proof, see Remark 15.31 Similarly to the wave equation, also here we find 
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that the rate of weak convergence is twice that of the strong convergence [17j under 
essentially the same assumptions. 

The literature on weak convergence for parabolic equations is richer. We have 
already mentioned 10 that proves results for fully discrete schemes of the linear 
stochastic heat equation. In |14j spatial, finite element schemes are considered for 
the same equation. Such schemes are also studied for the linear Cahn-Hilliard-Cook 
equation as well as the linear heat equation in [15J . Semidiscrete temporal schemes 
are investigated in [12] for the linear heat equation and in [S] for the nonlinear 
heat equation in one dimension. The techniques of the latter paper are extended to 
spatially semidiscrete schemes in multiple dimensions in pQ. A recent paper, |20) . 
successfully uses the methods of [10 to study a linear parabolic SPDE driven by 
impulsive noise instead of a Wiener noise. This indicates the possibility of extending 
the results of the present paper to larger classes of noise. 

2. Preliminaries 

Here we collect some background material from infinite-dimensional stochastic 
analysis and stochastic PDEs. We use the semigroup approach of DaPrato and 
Zabczyk and we refer to the monograph [5] for details and proofs. 

Let U and % be real separable Hilbert spaces; we often denote both their norms 
and scalar products by || • || and (•, ■) when the meaning is clear from the context. 
We denote the space of bounded linear operators from U to % by B(U, 7i) and the 
p-.th. Schatten class of operators from U to 'K by jC p (U, %). They are Banach spaces 
for all integers p > 1 and we will denote their norms by || • ||£ (u.t-C)- The operators 
in £i(W,H) are also refered to as trace class operators and operators in £ 2 (^,'H) 
as Hilbert-Schmidt operators. The space £2(M,'H) is a Hilbert space with inner 
product denoted (•, ■) c 2 {u,u)- When the underlying Hilbert spaces are understood 
from the context we will write || • || Tr = || • \\d(UM), II • IIhs = II ■ \\c 2 (UM) and 
(") - )hs = (•, •) c 2 {u.u) i n order to - we hope - increase the readability of the paper. 

In case U = hi we write B{U) = B{U,U) and C P {U) = L P {U,U) for short. If 
T 6 C\(U) and {etj^Lx is an orthonormal basis of U, then the trace of T, 

oo 

Tr ( T ) :=X^ Tefe ' efe ^' 
fc=l 

is a well defined number, independent of the choice of orthonormal basis. Below we 
state a number of properties of Schatten class operators. For proofs and definitions 
we refer to, for example, [6j Appendix C], [18] and [26] . 
If T e C P {U,H), then its adjoint T* e C P {U,U) and 

(2-1) \\T\\ Cp (UM) = \\T*\\c p (u,u)- 

If U = H and p = 1 , then also 

(2.2) Tr(T) = Tr(T*) 
and 

(2.3) |Tr(T)|<||T|| Tr . 

Further, if T is selfadjoint and positive semidefinite, then Tr(T) > and (|2.3j) holds 
with equality. 
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If Hi, Hi, and H are separable Hilbert spaces and T G C P (H 2 ,'H) and if Si G 
B(JAi,U 2 ) and S 2 G B(H,Ui), then 

, . 1 1 TSi I k P («!,«) < \\T\\c p (u 2 ,U)\\Si\\b(Ui,u 2 )' 

II^TIU,,^,^!) < l|r||£ p ( W2: «)||5 , 2 ||e('H,w 1 )- 
If S G £("H,W) and T G Ci{U,U), then we also have 
(2.5) Tr(TS) = Tr(ST). 

Moreover, iiT:U^U and T*T G Ci(U), then T G C 2 {U,U), TT* G £i(%) and 

HS 



T*T\\ Tl = Tr(T*T) = ||T||^ S = ||T*" 2 



( • ) =Tr(TT*) = ||TT*|| Tr . 

Finally, we note that if T G C 2 {U,U) and S* G C 2 {U,U), then TS G £i(H) and 

(2.7) ||TS|| Tr < ||T||hs||S||hs = (Tr(TT*)Tr(SS*)) 1/2 - 

To be able to compare various assumptions on the regularity of the noise, where 
the regularity usually is measured in the trace or Hilbert-Schmidt norms, we cite 
Theorem 2.1 in [15]. 

Theorem 2.1. Assume that Q G 23(%) is selfadjoint, positive semidefinite and that 
A is a densely defined, unbounded, selfadjoint, positive definite, linear operator in 
% with an orthonormal basis of eigenvectors. Then, for s G K, a > 0, we have 

\\AiQi\\i s < \\A S Q\\ TT < \\A s+a Q\\ Bm \\A-«\\ Tr , 

(2.8) \\AiQ^\\l s <\\A^QA^\\ Tl .. 

Furthermore, if A and Q have a common basis of eigenvectors, in particular, if 
Q = I , then 

\\AiQi\\%s = \\A s Q\\ Tl = \\A s +^QA-^\\ Tl . 

Let (fi, T, P) be a probability space and L p (yi,Ti) denote the space of random 
variables X : (il, J 7 ) — > (%, Bor(%)), where Bor(%) denotes the Borel er-algebra of 
the separable Hilbert space H, such that 

I!^IIl p (^)=E(||X|| p w )- / ||JC(o;)||? £ dP( W )<oo. 

By "strong convergence" we mean norm convergence in ^(^1,71). 

If {F(t)} t £[QT] is a family of (deterministic) bounded linear operators from a 
Hilbert space U to another Ti, then the Ito integral (also called the Wiener integral 
as the integrand is deterministic) 

/ F(t)dW(t) 
Jo 

with respect to a U- valued Q- Wiener process is well defined if 

(2.9) f \\F(t)Q l/2 \& s dt= f Tr(F(t)QF*(t))dt<oo. 
Jo Jo 

If (|2.9p holds then we have Ito's isometry 



(2.10) j T F(t)dW(t) ' = f \\F{t)Q^\\l S <tt- 
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The functionate G in (|1.4j) are called test functions and throughout this paper 
we will assume that they are mappings from H to R with bounded and continuous 
first and second Frechet derivatives. That is, they belong to the space 

Cl = Cl{HM) = {G g C 2 (H,R) : ||G|| c = (KiK) < 00} , 

where 

IIGUcg^R) := sup \\G'{x)\\ H + sup \\G"(x)\\ B(n) . 

The derivatives G'(x) and G"(x) are identified with elements in % and B("H), 
respectively, by the Riesz representation theorem. Note that || • ||c 2 ( , h,r) i s om y a 
seminorm and that we do not assume that G itself is bounded. 

3. AN ERROR REPRESENTATION FORMULA 

Following |10j and |15| , we start by developing a representation of the weak error 
e(T) in (|1 .4[) in the abstract setting. To do this we use the process {Y(t)}t>o in 
(|1.5p . being the unique weak solution of the drift-free differential equation (|1.8p 
with the important property that Y(T) = X(T). We also introduce the auxiliary 
problem 

dZ(t) = E(T -t)BdW(t), te(r,T]; Z(r) = f, 

where £ G Li(fl, H) is an Jv-measurable random variable. Its unique weak solution 
is given by 

(3.1) Z(t,T,Z)=S + j E(T-s)BdW(s), t£[T,T]. 

We note that Z(t, 0, E(T)X ) = Y(t) and that Z(t,t,£) = £. For G G C£(H,R), 
we define the continuous function u: H x [0, T] — >■ R by 

(3.2) u(x,t) = E(G{Z(T,t,x))). 

It follows from (|3.1|) and (|3.2|) that the partial derivatives of u are given by 

«*(*,«) =E(G'(Z(T, *,*))), 

Ura (x,t)=E(G"(Z(r,t,a;))). 
Hence, 

sup \\u x (x,t)\\ < ||G|| C 2 (KR) , 
(x,t)eWx[o,T] 

sup ||uxx(a;,t)||B(w) < ||G|| C 2 (W)R) . 
(x,t)eHx{o,T\ 

It is known that u is a solution to Kolmogorov's equation 

(3.5) 

u t (x, t) + \ Tr (u xx (a;, *)S(T - t)BQB*E{T - t)*) = 0, (x, t) e H x [0, T), 
u(x,T) — G(x), xeH, 

see, for example, Lemma 6.1.1], where we note that G G G^(H,R) is enough 
for existence and G G UC^(H, R) is only needed for uniqueness as the proofs of [7J 
Theorems 3.2.3 and 3.2.7] show and the global boundedness of G is not needed by 
Remark 3.2.1]. Finally, the partial derivatives of u in (|3.5j) are continuous on 
[0,T)xH. 

Our key result is the following representation formula for the weak error. 
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Theorem 3.1. Assume that ()1.2[) and (|l.ll|l hold and let {X(t)} t £[o,T] be the 
unique mild solution (|1.3j) of (jl.ip and that X(T) can be represented as X(T) = 
Y(T), where Y is given by (|1.10j) . 

If G G C^("H,M), then the weak error e(T) in (|1 .4[) /ias i/ie representation 

e (T) = E / ( U:c (y(o) + s(y(o)-y(o)),o),F(o)-y(o))ds 

(3-6) 7 ° ' 

+ iEy Tr(« xx (F(i), *)©(*)) di, 

(3.7) 0(t) = (-E(T - i)B + E{T - t)B)Q(E{T - t)B - E{T - t)B)* , 
or 

(3.8) 0{t) = (E(T - t)B - E(T - t)B)Q(E(T - t)B + E(T - t)B)* . 
Proof. As in [6l Theorem 9.8], since £ is J^-measurable, we have that 



(3.9) «(£,*) = E(G(Z(r,i,0) 

Thus, by the law of double expectation, 

(3.10) E(«(e,t)) = e(e(g(Z(T,£,0)|^)) = E(G(Z(T,t,£)) 
Therefore, taking also into account that X(T) = Y(T), it follows that 

e(g(x(t))) = e(g(F(T))) =e(g(Z(T,0,Y(0))) =e(u(Y(0),0) 
and, since X(T) — Y(T), we also have that 

e(g(X(T))\ =F,(g{Y{T))\ =F,(g(Z(T,T,Y(T))) =b(u(Y(T),T) 
Hence, 

e(T) = e(g(A(T)) - G(X(T))) - e(u(F(T),T) - u(F(0),0)) 
= E(u(Y(0), 0) - u(Y(0), 0)) + e(u(F(T), T) - U (Y(0), 0)) . 
For the first term we note that due to the differentiability of u we can write 

E( M (f(o),o)- u (r(o),o)) 

= E j (u x (Y(0) + s(Y(0) - F(0))),F(0) - F(0)) ds. 

For the second term, we use Ito's formula (see [U Theorem 2.1], where, in contrast 
to [6l Theorem 4.17], uniform continuity of the appearing derivatives on bounded 
subsets ofHx [0, T) is not assumed) for u(Y(t),t) on [0, T — e] and passing to the 
limit e — > 0+ using the continuity of u on T-L x [0, T] and the continuity of the paths 
of Y{t) on [0,T]. Thus, taking also Kolmogorov's equation (|3.5I) into account, we 
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get 



E(u(y(T),T)-u(y(0),0)) 

= e| {u t (Y(t),t) 
(3.11) + ± Tr (u xx (Y(t),t)[E(T - t)B]Q[E(T - i)S]*) } dt 

= ±E /" Tr ( W »(y(t),(){[£(T - t)B]Q[E(T - t)B}* 

- [E{T - t)B]Q[E(T - t)B]*}) dt. 



The operator u xx (t;,r) is bounded for every £ and r and both E(s)BQ[E(s)B]* 
and E(s)BQ[E(s)B]* are of trace class for almost every s by assumptions (II. 2\i and 
p. lip . Hence, the trace above is well defined for almost every t since by (|2.4I) with 
P=l, 

|K x (£,r)[£( S )B]Q[£( S )B]*||Tr < IK.(^)lls(«) ||[£( fl )S]Q[S( fl )S]*||iv 

= ||«»xK,r)||B( W ) Tk([£(«)B]Q[£(«)B]*), 

where the last step is (|2.3|) with equality, which holds since [-E'(s)B](5[-E l (s)B]* is self- 
adjoint and positive semidefinite. The same computation can be done with the term 
involving [E(s)B]Q[E(s)B]*. Further, the operator u xx (£,r)[E(s)B]Q[E(s)B]* is 
also of trace class for almost every s, since, by (|2.1j) . (|2.4j) . and (|2.7|l . 



||7/ xa; (£,r)[i?( 5 )i?]Q[i?( S )i?]*|| Tr 

< IK*(£,r)|| B(w) ||[S( S )i?]Q[S( s )B]*|| Tr 

< ||«xx(e,r)||B(«) |P(s)S]Q 1/2 ||hs HQ^TO^'IIhs 
- IK*(£,r)]] S(W ) II^^SlQ^n^ ||^( s)j B]QV2|| Hg 

= ll^(e,r)|| K( « ) (Tr([£;( S )S]Q[£;( S )i3]*)Tr([ J B( S )B]Q[^( S )S]*)) 1/2 . 

Therefore we may rewrite the operator in the trace in (|3.1ip by adding and sub- 
tracting u xx (£,r)[E(s)B]Q[E(s)B]* to get 

u xx (tr){[E( s )B]Q[E(s)B}* - [E(s)B]Q[E(s)B]*} 
= u xx (t r)[E(s)B - E(s)B]Q[E(s)B]* 

+ u xx (t r)[E{s)B]Q[E(s)B - E{s)B}* 
= :Oi+0 2 . 
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Further, using ()2.2[) . (|2.5j) . and that Q and u xx (£,r) are selfadjoint, we obtain 
Tr(Ox + 2 ) = Tir(Oi) + Tr(0 2 ) = Tr(Oi) + Tr(C>*) 

= Tr(Oi) + Tr([E(s)B - E(s)B}Q[E(s)B]*u xx (S, r)) 
= TV(0 1 ) + TV(u xx (e,r)[ J B(s)S- J B(s) J B]Q[ J B(s) J B]*) 

(3.12) = Tr ( Ubx (£, r)[B(a)B - £(s)B]Q[£(s)£ + £( a )B]« 

= Tr ([E(s)B + E( S )B]Q[E{s)B - E{s)B]*u xx ^, r) \ 

(3.13) = Tr (u xx (£, r)[E(s)B + E{s)B]Q[E(s)B - E(s)B]* 
Finally, by inserting (J3.12I) or (j3.13[) into (I3.1ip we finish the proof. □ 

4. Application to the wave equation 

4.1. A general error formula. In this section we apply the general result from 
Section [3] to the stochastic wave equation (|1.12p . As mentioned in the Introduction 
the wave equation may be re- written in the form (jl.ljl . We start out by making this 
statement precise. At the same time we introduce a framework for measuring the 
regularity of the solutions and to perform a careful error analysis. To this aim we 
let V denote a convex open bounded domain in M. d with polygonal boundary dT>, 
and equip L 2 (T>) with the usual norm || • \\l 2 {v) an d inner product (•, ■)l 2 (d) an d 
let A = -A be the Laplace operator with D(A) = H 2 (V) n Hq(D). To measure 
regularity we introduce a scale of Hilbert spaces. Let {(\j, 4>j)}j^i be eigenpairs 
of A with a nondecreasing sequence of positive eigenvalues Xj and corresponding 
orthonormal eigenvectors <f>j. For a£l we endow D(A a ' 2 ) with inner product 

oo 

(v,w) Aa := {K a/2 v,K a/2 w) L ^ (v) =^A J Q («,^) L2(P) ( W , ( /) J ) L2(I , ) , 

i=i 

and the corresponding norm ||v||?v a = (v,v)^ a . For a > the space H a now may 
be defined through 

(4.1) H a := {v e L 2 {V) : ||«|| Ao < oo}. 

If a < 0, then H a is taken to be the closure of Li{J5) with respect to || • ||^ Q . 
It is notable that with a > the space H~ a may be identified with the dual of 
H a and that H a C H? if a > (3. Further, it is known that H° = L 2 (V), H 1 = 
H$(T>), H 2 = H 2 (V)r\H^{V), see [24j Chapt. 3]. In addition we define the product 
space 

(4.2) U a := H a x H"- 1 , a e R, 
with inner product 

(V,w) n <* = {V\,Wx)jj a + (V2,V)2)h«-i> 

where v — [wi,W2] T and w — [wi,W2] T - The corresponding norms are 

\\v\\ 2 n » :=ll«i||^ + 11^11^. 
We take H to be the special case of (|4.2[) when a = with norm || • || = || • ||^o and 
inner product (•, •) = {-,-)l 2 (V) + (-A- -1 ' 2 ") A~ 1 ' 2 ')l 2 (v)- We now regard A as an 
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operator from H 1 to H 1 as (Ax)(y) = (Vx, ^7y)L 2 (D) and let A and B be defined 

by 



(4.3) 



A:-- 



-I 

A 



D 



Here _B is considered as an operator from H 1 to H and the domain of A is 



D{A) = Ixen-. Ax = 



-x 2 
Axi 



G U = H" x 27 



7- 1 }=^ 1 = ff 1 



xff u 



We have some freedom in defining U but a natural choice is £Y = H° = L2CD). 
Thus, we consider the process {W(t)} t >o to be a Q- Wiener process in L2CD) and 
thus Q to be bounded, selfadjoint and positive semidefinite on U = L^iT)). Note 
that H° C i? -1 , and therefore -B G BiU^H). It is well known that the operator 
—^4 is the generator of a strongly continuous semigroup E(t) — e~ tA on H, in fact, 
a unitary group, that can be written as 

(A A\ F(t\ - P - tA - \ C ^ *- 1/2 S(t) 

(4-4) E(t)-e -[_ A i/2 S{t) c{t) J, 

where C(t) = cos(tA 1 / 2 ) and 5(4) = sin(tA 1 / 2 ) are the so-called cosine and sine 
operators. 

With these definitions (jl.ll) becomes the stochastic wave equation (|1.12l) with 

solution [Xi(t), AT 2 (i)] G W if also the initial value X = [X , 1,^0,2] G H and 
dOJ holds. 

We are know ready to use the framework set forth in Section [31 starting by 
assuming as little as possible about the type of perturbations of X(i). To get a result 
that enables error analysis of the full solution X, as well as either of the coordinates 
X\ or X<i, we make three additional, rather weak assumption on the data of the 
problem. First, we take G to be the composition of a function g G C^(V, R), where 
V is a real separable Hilbert space, with an operator L G B(H, V), i.e., 

G(v)=g(Lv), v = [v 1 ,v 2 f eH. 

It is clear that G e Cg(M,I). Typically we have V = U and L = I, or V = H° 
and Lv — P x v = v\ (projection to the first coordinate), or possibly V = H^ 1 
and Lv — P 2 v = V2 (projection to the second coordinate). Second, since in all 
interesting cases the function Xq will be related to Xq we will assume that this 
relation can be described by an operator P G B(H) through 

Xq = PXq. 

In the cases we consider below, P will be the orthogonal projection onto a finite 
element subspace of % or the identity operator but it could also be an interpolation 
operator. The third assumption we make is that B = PB. This is unnecessarily 
restrictive in general but it suffices for the purposes of this paper and it also makes 
the presentation clearer. We begin with a technical result. 



Lemma 4.1. If E(t) is given by 
equivalent. 

(i) Tc{hr l / 2 Qhr 1 / 2 ) < 00. 
(ii) ||A- 1 /2qi/2|| hs<0O . 

(iii) HA-Vsqa- 1 /*!!^ < 00. 



then the following four statements are 
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(iv) Tr(/ T E(t)BQB*E(t)* dt) < oo for some, hence all, T > 0. 
If either of them holds, then 

(4.5) Tr(f E{t)BQB*E(t)*dt) ^TTt{A-^ 2 QA-^ 2 ). 

Proof. The operator A _1//2 QA -1 / 2 is selfadjoint and positive semidefinite on H°. 
Hence, if the trace is finite it equals the trace norm and the other way around. 
This implies that Q <^> (fm]) . We know assume that ([Iv]) holds. By monotone 
convergence, 

(4.6) Tr ( I E(t)BQB*E(t)* dt) = f Tr (E(t)BQB*E(t)*) dt 

and by flUJ) and since E(t)* = E(t)~ l we have 
Tr n (E{t)BQB*E(t)*) = Tij ia {Q 1 ' 2 B*E{t)*E{t)BQ 1 ' 2 ) = r Ii f{Q {Q l ' 2 B*BQ 1 ' 2 ) 

— \\BQ \\c 2 (HO,-H) = HI"' ^ J IU 2 (H°,H) = 11^ \\c 2 (H°,H- 1 ) 

= ||A- 1/2 Q 1/2 || 2 £2(i j 0) =Tr^ ) (Q 1 /2 A -iQi/2 ) =T r^(A- 1 /2Q A - 1 /2). 

Therefore, it follows that the integrand on the right-hand side of (|4.6[) is constant 
and is equal to Tr J j (A~ 1 ' /2 <3A _1 <'' 2 ), which implies (|4.5j) . Therefore © must be 
true. This argument is reversible so Q 4=> (JIvJ). But from this computation it is 
also evident that 

I|A- 1/2 Q 1/2 Hhs = TrtA-^QA- 1 / 2 ); 
that is, that <£> ©. □ 

Corollary 4.2. J/ | A' 3_1 / 2 QA _1 / 2 || Tr < oo for some (3 > 0, tften (EvJ) ZioZds and 
(jl.ip /ias a unique weak solution {X(t)}te[o T\- 

Proof. By pi5]) and (J2H) we have that 

|Tr(A- 1 / 2 gA- 1 /2)| = HA-^QA- 1 / 2 !!^ = HA-^-^QA- 1 / 2 !!^ 
< [[A-^Ha^jHA^-VaQA-Vaii^. 
The result now follows by Lemma [4.11 □ 



Since HA^gl^ < | A' 3 - 1 / 2 QA- 1 / 2 || Tr by ([H]), we may conclude that under 
the assumption in the previous corollary, we actually have mean-square regularity 
of order (3. This follows from the following theorem, which is quoted from [TBI 
Theorem 3.1]. 

Theorem 4.3. If ||A^Qi ||fj g < oo, X G L 2 (n,H p ) for some (3 >0 and A,B 
as in (|4.3p . i/ien i/ie mi/d solution (|1.3j) satisfies 

\\X(t)\\L 2 (n,nn < C(\\X \\ La{a , m +T 1 / 2 ||A^ i Q^||Hs). 

We will need the following result on the Holder continuity of E(t). It will put 
an ultimate limit on the convergence rate that one can achieve with respect to 
time-stepping. 

Lemma 4.4. If {E(t)}t>o is the semigroup in (|4.4[) . then 

\\(E(t) - E(s))x\\ <C\t-s\ a \\x\\ n c, xGH a , t,s>0, ae [0,1]. 
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Proof. The operator E(t) is bounded on % so the statement is true for a = 0. Let 
a = 1 and x — [xi,x 2 ] T G H 1 . Then 

\\(E(t) - E(s))x\\ 2 = ||(C(t) - C(*))a:i + (5(t) - S(s))A-^x a \\%o 

+ ||(S(») - S(t))A^ a xi + (C(i) - C(«))a !a ||^_ 1 =: A 1 + A 2 . 

By the triangle inequality and the definition of the norm on H , we have for the 
last term that 

A 2 < 2||A- 1 / 2 (5(t) - S(s))A 1 ^ 2 x 1 \\ 2 {0 + 2\\A- 1 / 2 (C(t) - C{a))x a \\%o 

= 2\\(S(t) S(s))xi\\% + 2\\(C(t) C{ S ))K-^x 2 \\\ a . 

Since x 2 G H° it follows that A~ 1 / 2 x 2 G H 1 . Hence we only need to investigate the 
Holder continuity of C and S as functions from [0,T] to B{H l , H°). To this aim 
we note that for real y the inequality 

(4.7) \sm.{ty)-sm{sy)\<\t-s\\y\ 

holds. It follows that for £ e H 1 we have ||(5(t)-5(s))^||^ < |*-s|||f H^. Indeed, 

OO 

\\(S(t) s(aM\\%, = E« 5 W ~ 5 ( s ^ &>*> 

3=1 



J2( S m(tXy 2 )-8m(sXy 2 )) 2 ^4>j)] 



) 2 

3=1 



3=1 

The inequality (I4.7|) holds also with sin replaced by cos. Thus the statement of the 
lemma holds also for a — 1. The intermediate case follows by interpolation. □ 

We are now ready to prove a weak error bound for perturbations of the stochastic 
wave equation. 

Theorem 4.5. Assume that {A(£)} te [ ,T] *s the mild solution p.3[) of the stochastic 
wave equation ifTTj) ««#» ||A' 9 ~ 1 / 2 <9A~ 1 ''' 2 || Tr < oo and X G Li(£l,% 2 ' 9 ) /or some 
/3 > 0. Assume further that X(T) can be represented as X(T) = Y(T), where Y(t) 
is given by (fTTUl) with X = PX , P G B(H) and B = PB and such that ([TTT]) 
holds. Let g 6 Cg(V,R) and L G B(H,V). Define 

(4.8) ^i:= sup \\E(t)P\\ n , 

te[o,T] 

(4.9) i^ 2 := ||A' 9 - 1 / 2 QA- 1 / 2 || Tl .. 

Then there is C = C(T, \\L\\ B{Hy) , ||^'o|U 1 (n,'H 2 ' 3 )) \\g\\c^ K i^ K 2) such that 

(4.10) v(g(LX(T))-g(LX(T))) <C sup \\L(E(t)P - E{t))\\ B{w , >vy 

v ' te[o,T] 

We want to emphasize that this theorem reduces the problem of proving weak 
error estimates for the stochastic wave equation to proving deterministic error es- 
timates for the approximation of the semigroup or, to be precise, to find a bound 
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sup \\L(E(t)P-E(t))\\ B(n2?tVy 
*e[o,T] 



for 



Proof. First note that (|4.9j) implies (|1.2|) by Corollary 14.21 and hence we may use 
TheorernOwith G(X) = g{LX), To this aim we note that G'(X) = L*g'{LX) and 
G"(X) = L*g"{LX)L. The terms in (|3.6[) will be estimated in order of appearance 
with O as in (JSU). For the first term, by (O, (O and the fact that both F(0) 
and Y(0) are J-Q-measurable, we have 



E 



u x (Y(0) + s(Y(0) -Y(0)),0),Y(0) -Y(0))ds 



(g'{LZ(T, 0, Y(0) + s(Y(0) - Y(0)))) \t ) , L(Y(0) - Y(0)) 



< S u P \\g'(x)\\E(\\L(Y(0)-Y(0))\\) 

= sup \\g'(x)\\ E(\\L(E(T)X - E(T)X )\\) 
xev 

< \\g\\ cl \\L(E(T)P - EiTm^^^Xoh^^y 

For the second term we note that with 

0~(t) = E(T - t)B - E(T - t)B = (E(T - t)P - E(T - t))B, 
+ {t) = E(T - t)B + E(T - t)B = (E(T - t)P + E(T - t))B, 

and by (|3.3[) . (J3.9I) . we may write 
E 



d.s 



Tr [u xx (Y(t),t)0 + (t)QO-(t)*j di 

= E J Tr(v{L*g"(LZ(T,t,Y{t)))L\Tt)0 + (t)QO-{t)*\&t. 

Using (|2.5p . (I2.3P and (|2.4[) we bound the integrand above as follows: 
Tr (e(L* g" '(LZ(T,t,Y(t)))L\^ t )0 + (t)QO~ ■(«)*) 

= Tr(E(. 9 "(LZ(T,t,f(t)))|j' t )iO + (t)gO-(i)* J L*) 
< sup ||.g"(x)|| B(v) ||LO+(i)QO-(t)*L*|| T ,. 

Here we have sup xeV ||s"(a;)||;B(v) < IMIc 2 an d by (|2.1|) and (|2.4|l . 

||iO + (t)QO-(t)*L*|| Tr = [[iO+(t)A 1 / 2 A- 1 /aQO-(t)*L*|| lk 

< HLO+WA^II^o >v) HA-Vago-^U*!^ 
= || J LO+(t)A 1 / 2 H S(A o >v) HLO-CtJQA-Vail^ 

= HLO+WA^Hb^o , v) HXO-WAVa-zSA^-VaQA-Vai^ 

< ||iO + (i)AV2|| B(i . 0]V) [[iO-^AVa-zs \\ B{ Ho y) \\A^ 1/2 QA- 1/2 H. 
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The first factor can be estimated as 

WLO+WA^Wb&m = \\L(E(T - t)P + E(T - t))BA l / 2 \\ Bi HO, v) 

< ||L|| B(?iiV) (||£(r - t)P\\ B(n , n) + \\E(T - t)\\ B{nM) ) \\BK x t*\\ B ^ w 

< ||i|| B(w .v)(^i + l), 

because \\E{s)\\bchm) = 1 = W^^^Wbih HY Similarly, the middle factor may be 
bounded by 

\\L(E(T - t)P - E(T - t)) \\ B(H „ y) \\BA 1/2 - I3 \\ b{ ho^ } , 
where \\BA 1 / 2 ~^\\ B ,j i0 u - lfi - ) = 1. The third term is K2. Thus, we conclude that 

E / Tr (u xx (Y(t),t)0-(t)QO + (t)*) dt 

<C^ \\L(E(T - t)P - E(T - tj) || B(wa/J|V) dt 

<CT sup \\L(E(t)P - E(t))\\ 

te[o,T\ y ' ' 

and the proof is complete. □ 

Remark 4.6. We briefly comment on the abstract condition \\A i3 ~ 1 / 2 QA~ 1 / 2 \\tt < 
00 of Theorem 14.51 If Q = I, then we must have Tr(A^ _1 ) < 00 and hence taking 
the asymptotics of the eigenvalues of A into account we conclude that §03— 1) < — 1; 
that is, (3 < 1 — |. Hence d = 1 and /3 < h. In general, using (|2.4j) . we have that 

WAP-WQA-WWn < Tt^-Va-JHA'QA-*!^. 

Thus, if |(jS - I - s) < -1; that is, /3 < s + \ - f , then HA^-VaQA- 1 ^]]^ < oo 
provided that A S QA~2 is a bounded linear operator on U. 

4.2. Weak convergence of temporally semidiscrete schemes. We begin by 
applying Theorem l4.5l to semidiscrete approximation schemes where the discretiza- 
tion is with respect to time. We will use results on so-called /-stable rational 
approximations considered in [3]. An /-stable rational approximation of order p 
is a rational function R such that (|1 .13|) holds. A class of such functions are con- 
structed in [2] and analyzed further in connection to oscillation equations in |21) . 
It contains the implicit Euler method (p = 1) and, which is important since it 
preserves energy for the wave equation, the Crank-Nicolson method (p = 2). 

For these functions the operators E). — R(kA), k > 0, are well defined on H and 
they are contractions, and hence stable, as —A generates a unitary group. Here 
k = T/N, N G N, is the time step. We can approximate the solution of (|1.1|) on 
the uniform grid tj — jk, j = 0, . . . , N, by the solution of the difference equation 

(4.11) X{ = EkiXf 1 + BAWl), j = l,...,N- X° = X , 

given by 

X™ = E%X + J2 E%- j+1 BAW j , n > 1 , 
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where AW j = W(tj) — W(tj-i). We want to define a process {Yfc(t)}te[o,r] of the 
form (|l.fO|) that is as close as possible to (|1.5p and such that Xj? = Y k (T) . To this 
aim we first define a new discrete process 

n 

(4-12) Y£ = E*~ n X% = EJ?X + J2 E^- j+1 BAW j . 

3=1 

Clearly Y^ = X^ . In order to make a piecewise constant time interpolation of 
(|4.12[) we introduce the time intervals Ij — [tj-i,tj) for j = 1, . . . , N and In+i — 
{i./v} = {T}. With x being the indicator function we then write 

JV+l 

E k {T-t)=Y J E^ +1 Xlj {t). 

3=1 

It may easily be checked that this corresponds to writing 

N 

E k (t)=j2 E ixt j (t) 

3=0 

with Jo = {to} = {0} and Ij = (tj—i,tj] for j = 1, . . . , N. We finally define 

(4.13) Y k (t):=E k (T)X + f E k (T - s)B dW(s). 

Jo 

The process Y k has the desired properties and, in addition, Y k (t n ) = Y k n . To apply 
Theorem 14.51 it remains to show that p. lip holds with E(t) — E k it) and P = I 
(hence B = B). This is indeed the case as soon as we are guaranteed a weak 
solution of (jl.ip . as stated in the following lemma. 

Lemma 4.7. If E{t) = E k (t) and B = B, then (fOJ) implies IjTTTTt . 

Proof. If dO} holds then the trace of A" 1/2 QA- 1 / 2 is finite by Lemma |4~T] Thus, 
for alii > 0, using (j2~3)l and (f2~4l) . it follows that 



Tv(E k (t)BQB*E* k (t)) < \\E k (t)\\ 2 B{n) \\BQB*\\ Tl 

< \\BQB*\\ Tt = Tr(BQB*) = Tr^-^QA- 1 / 2 ) < oo, 
where the last equality is shown in the proof of Lemma 14.11 as 

Tr(BQB*) = Tr(Q 1/2 B*BQ 1/2 ) 
by (|2.6p . The statement of the lemma follows by the monotone convergence theorem 



again as in the proof of Lemma 14.11 □ 

In order to make use of Theorem 14.51 we need a bound on E(t) — E k (t). The 
results in [3] are concerned with the difference at the grid points. With our notation 
their conclusion reads that, for x G H p+1 , 

(4.14) \\(E(t n ) - E%)x\\ < Ct n kP\\x\\ HP+1 . 

As already mentioned, the conditions in (11.131) ensures that the operator E% is a 
contraction on H for any n > 0, so (|4.14p can be extended to fractional order by 
interpolation, i.e., 

(4.15) \\(E(t n )-E%)x\\<Ct n k a T&\\x\\ n ., ae[0,p+l]. 
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For our purposes, it is not enough to consider only the grid points, but fortunately 
a global error estimate follows easily. 

Lemma 4.8. For the operators E k (t) and E(t) defined above, we have that 

sup \\E k (t) - E(jt)\\tn H c tH) < CCOfc™^' 1 ), k > 0, 
te[o,T] 

where p is a nonnegative integer as in (|1.13|) and a > 0. 

Proof. The statement of the lemma follows from (|4.15|) and Lemma 14.41 Indeed, 
for t G Ij, we have 

E k (t) - E(t) = (E k (tj) - E(tj)) + (E(t,) - E(t)) = (E{ - E{t,)) + (E(t,) - E(t)). 

Hence, with I = [0, T] and 1 = {0, 1, ... , N}, 

sup\\E k (t)-E(t)\\ B{naM) 
tei 

< sup (\\El ~ EMWbch-.h) + ™p \\E(tj) - E{t)\\ B(Ha . H) ) 

< C(T)(fc min(Q 5*T' p) + fe min(a,l)) < ^^minfa^l) ^ fc < L 

Finally, for fc > 1, the statement follows by stability. □ 

We are now ready to prove a bound for the weak error of the pure time- 
discretization via (|4.11[) of the stochastic wave equation. 

Theorem 4.9. Assume that HA^-^QA- 1 / 2 ^ < oo and X G Li(fi,W 2/3 ) for 
some (5 > and G G C 2 ('H,M.). Then the weak error of the rational approximation 
algorithm (|4.11|) of the stochastic wave equation described above is bounded by 

(4.16) \E(G(X?)-G(X(T)))\ < Ck min{ ^^' x \ fc > 0. 

Proof. We may use Theorem S3] with E(t) = £? fc (t), P = I, V = U and L = I, 
because ||A /3 ~ 1 / 2 QA~ 1 / 2 || T r < oo implies (jl.lll) by Corollary I47Z1 and Lemma 14771 
Thus, our claim follows by applying Lemma l4~8l with a — 2/3 to (|4.10[) . D 

4.3. Weak convergence of fully discrete schemes. In this section we will 
present an error estimate for a fully discrete scheme. We will borrow the set- 
ting from [2], where estimates for the deterministic wave equation are proved. The 
spatial discretization is performed by a standard continuous finite element method 
and the time discretization, as above, by /-stable rational approximations of the 
exponential function. We briefly describe this method and state the error estimates 
from [2]. 

We assume that T> is a convex polygonal domain and we let {SX}o</i<ij r — 2, 3, 
be a standard family of finite element function spaces consisting of continuous piece- 
wise polynomials of degree r — 1 with respect to a regular family of triangulations 
of V. Moreover, we define S r hfi = {»6SJ: v\ dv = 0}, so that S r hQ C H 1 . With 
HP denoting the standard Sobolev space we then have the error estimate 

(4.17) \\R h v-v\\<Ch \\v\\ H P, vefrnH?, /?e[l,r], 
where the Ritz projection Rh : H 1 — >• S^ is defined by 

(VR h v, V X ) - (V», V X ), to e ff 1 , * e S r hfi . 
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Further, we define the discrete Laplacian Ah : S^ Q — > S^ by 

(A h 7 7 ,x) = (V? 7 ,Vx), Vj),x£^ . 
The homogeneous spatially semidiscrete wave equation is to find 
u h (t) := [u h>1 (t),u h , 2 (t)] T G ST x ST 



such that 
(4.18) 



Uh.2_ 


+ 


' -I" 

A h _ 




W/i,l(*) 

U h ,2(t)_ 


= 


"o' 





, *>0; 



'ttfc,i(0)' 

«m(o) 



i\lU0,i 

Ph,2U ,2 



Here P/^i : iJ° — >■ 5£ and Ph,i ■ H l — > S^ are the orthogonal projectors defined 
by (P h ,if, X) = (/, x), V X G SI , for / e H° if i = 1 and / e H~ l if i = 2. 



It is well known that A/j has eigenpairs {(<f>h,j, ^-h,j)}j=i, where {A/ t j},-_' l 1 is a 
positive, nondecreasing sequence and {0/j,j},-Ji an 7J°-orthonormal basis of S^ lQ . 
If we write 

" -I] 



A h := 



A, 







and if Ph = [Ph,i> Pfi,2] T and uq := [uo,i, u o,2] T , then (I4.18P may be written 

(4.19) u h + A h u h = 0, t > 0; u fc (0) = P/^o- 

The operator — A/, is the infinitesimal generator of a strongly continuous semigroup 
Eh(t) and the solution of ()4.19|) is given by 

u h (t) = E h (t)P h u . 

Similarly to (|4.4[) the operator Eh (t) has a representation in terms of sine and cosine 
operators; i.e., 



E h {t) 



^h 1/2 s h (t) 



c h (t) 

-A h /2 S h (t) C h (t) 



with S h (t) = sm(tAl /2 ) and C h {t) = cos(tA* /2 ). 

The time discretization, as in the previous subsection, is performed by /-stable 
rational single step schemes; i.e., schemes where the rational function R fulfills 
(|1.13p for some positive integer p. The fully discrete problem on the same uniform 
grid as in Subsection 14.21 then reads 



(4.20) 



vl k = R(kAh)v n -\ n = h 



,N: 



J h,k 



PhUo- 



We will henceforth write Eh,k 
written as 

(4.21) 



R(kAh) and the solution of (|4.20|) may then be 



J h,k 



= El k P h u . 



The error estimate proved in [2] is as follows. It provides only a bound for the first 
component in u, which we express by means of a projector P . 

Theorem 4.10. If P 1 : H -> H° is defined as P 1 x = x\ for x = [x 1 ,x 2 ] T G U, 
then 

\\P\El k P h - E(t n ))u \\ A o < C(t n )(h r \\u \\ nr+1 + kP\\u \\ nP+1 ), t n = nk> 0. 

Using the stability of E(t) and EJ^ k and a standard interpolation argument, this 
results in the following bound on the error operator. 



WEAK CONVERGENCE OF FULLY DISCRETE APPROXIMATIONS 19 

Corollary 4.11. Under the assumptions of Theorem \4-.10\ we have, for /3 > 0, 

\\P\El k P h -E{t n ))\\ B{w<k0) < C{t n ){h^^^r) +k mM fil ^, P )^ ^ = nk > 

We return to the stochastic wave equation whose fully discrete version now reads, 

.,N; X°=P h X . 



with Bh := 


= P h B = 


[0,P h ,2] T , 






(4.22) 


x l,k 


= E hA X h7k 


+ B h AW j ), j = l,. 


The solution is given by 






(4.23) 




X h,k = E h l k P h X 


n 



B h AW J . 

3 = 1 

As in the previous section we multiply by E h 7 n and arrive at the drift free version 

a 

Yhm = EfoHXo + Y, <u 3+1 BuWi 
and with piecewice constant interpolation 

(4.24) Y htk (t) = E h . k (T)P h X Q + f E Kk {T - s)B h dW(s) 

Jo 

in exact analogy with the temporally semidiscrete case in (14.131) . 

Next we bound the weak error for fully discrete schemes given by (|4.22[) . We 

only prove a result for the first component in X . 

Theorem 4.12. Assume that ||A' 9 - 1 / 2 QA- 1 /2|| li , < oo and X G L^U 2 ?) for 
some /3 > 0. If X£ k is given by (|4.20[) and X(t) is the weak solution p. 31) of (II. ip 
with A, B as in (14.31) . then for g G Cg(H°,R), we have 

\E(g(x£ kA ) -g(X x (T)))\ < C{T)(h min{2 ' 3 ^TT> r) + fc"M 2 /»5$i>i)). 

Proof. The function in (j4~2^)) is clearly of the form ([HD|| with Y h}k (T) = X^ k and 
we have already seen that ||A^ -1 / 2 QA -1 / 2 ||Tr < oo implies (jl.21) . Furthermore, 

Tr(E htk (t)B h QB* h E h>k (t)*) < Tv(B h QB h ) < oo, 

as BhQBh is a bounded operator with finite-dimensional range and hence it is 
of trace class. Therefore also (jl.lip holds and Theorem 14.51 can be applied with 
V = H°, L = P 1 (as defined in Theorem glO]) , E(t) = E hik (t), B = B h and 
P = Ph. From Corollary 14.111 and Lemma [4.41 as in the proof of Lemma [4. 81 it 
follows that 

sup \\P\E htk {t)P h - E{t))\\ B{n2 , k0) < C{T) (h min ^^r^ + k min ^^^ 

t€[0,T] ' V 

Finally the statement of the theorem follows from inserting this into (|4.10[) . D 

4.4. Strong convergence of fully discrete schemes. It is a general phenome- 
non that the order of weak convergence is twice the strong order under the same 
regularity of the noise. This essentially turns out to be the case also for the sto- 
chastic wave equation discretized by the method described in the previous section. 
As we are not aware of any results on strong convergence of a fully discrete ap- 
proximation of the stochastic wave equation using finite elements in the spatial 
domain, we give a short derivation of the strong order. We remark that the strong 
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convergence is studied for a spatially semidiscrete finite element method in [16j , for 
a fully discrete leap-frog scheme in one spatial dimension in [25] , and for a spatially 
semidiscrete scheme in one dimension in [23 . 

First we form the strong error by taking the difference of (|4.23|) and (11.31) . pro- 
jecting onto the first component, and taking norms: 

E(||P 1 (^ fc - X(T))\\% ) < Cv(\\P\E» k P h - E(T))X \\%) 



C*E 



P 1 f (E h , k (T-s)P h -E{T-s))BdW(s) 
Jo 



H" 



h+h- 



If X eL 2 (fi,H /3 ), then 

h < C\\P\E h . k P h - E(T))\\l {H ^ Ao) E(\\X o f H0 ) 

by Corollary 14. Ill For J 2 we use Ito's isometry (12.10[) to get 



J 2 = e(| J P 1 {E Kk {T- S )P h -E{T-s))BdW(s) U2 



H'> 



= f \\P 1 (E Kk (T-s)P h ~E(T- S ))BQ 1 / 2 \\ 2 m d S 
Jo 

= f \\P l {E h , k (T - s)P h - E{T - S ))Bk^ K^Q^Wls&s 
Jo 

< £ \\P\E Kk (T - s)P h - E(T - S ))||^^ 0) ||A^Q 1 /2||2 is ds 
<T sup {\\P\E htk (t)P h - E{t))f k0) )\\l\^Q l/2 \\l s 



i-f 



where the first inequality follows from the fact that ||-BA 2 H^^o fjp\ = 1 combined 
with (|2.4[) , and the last inequality from Corollary 14.111 and Lemma 14.41 as in the 
proof of Lemma 14.81 Combining the bounds for I\ and I2 and taking square roots, 
we have shown the following result. 



Theorem 4.13. Let Wk—Q 1 /' 2 ^ < 00 and X Q e L 2 {fl,'H 13 ) for some ft > 0. 
Then the strong error of the approximation Xj^ k x — P l X^ k of the displacement 
Xi(T) = P 1 X(T) in the stochastic wave equation is bounded by 

lK fe)1 -Xi(T)|| i2(M0) < C(T)(h™V^ + fc"*^5*r. 1 )). 

The regularity assumption on Q in Theorem 14.121 implies the assumption in 
Theorem l4~T3l see Theorem |2~T1 with s = (5 - 1 in (l2~8|) . Thus the claim that the 
weak rate is essentially twice the strong rate is justified (if (3 is not too large) by 
comparing Theorems 14.121 and 14.131 Note also that the mean-square regularity is 
of order [3 according to Theorem 14.31 
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5. Application to parabolic equations 

Here, we give a detailed weak error analysis of a fully discrete scheme for the 
linearized Cahn-Hilliard-Cook (CHC) equation and also comment on the linear 
stochastic heat equation. 

The linearized CHC equation, see [17], is 

dX + A 2 A dt = dW, t> 0; X(0) = X , 

where now A = — A is the Laplacian together with homogeneous Neumann bound- 
ary conditions. To write the CHC equation in the form (|l.l|l we therefore set 
U = {/ £ L 2 (V) : {/, 1) - 0}, A - A 2 with D(A) = {/ <= H 2 (V) n U : §£ = 0}, 
where T> is a convex polygonal domain, and we take U = % and B = I. We further 
define the the spaces ff a = £>(A a ) in analogy with Section [4j Thus, H = H a , 
D(A) = H 4 and —A is known to be the infinitesimal generator of the analytic 
semigroup E{t) = e~ tA = e~ tA on TL. 

We recall the finite element spaces S^ (without boundary conditions) of order 
r = 2,3 from Subsection 14.31 and set S^ = {v E S^ : (v, 1} = 0}. We now define 
the discrete Laplacian A/j : S^ — > SJ^ and the Ritz projector Rh : H 1 — > 5£ in the 
analogous way and we have an error bound of the same form as in (J4.17I) . We set 
Ah = A 2 and note that — Ah is the generator of an analytic semigroup Eh(t) on 
Sh- We consider only the backward Euler time-stepping and therefore introduce 
Eh.k = (1 + kAh)^ 1 and define Eh,k(t) in an analogous fashion to the case of the 
wave equation, see (|4.13|) . 

We need error bounds for the approximation of the semigroup. We claim that, 
for all v £ H, 

(5.1) \\(El k P h - E(t n ))v\\ < C(h a + k a / 4 )t- a / 4 \\v\\, t n = kn, ae [0,r], 

where Ph'.'H—} S^ denotes the L2(^)-orthogonal projection to S^. To see this we 
write 

El k P h v - E{t n )v = (El k P h v - E h {t n )P h v) + (E h {t n )P h v - E{t n )v). 

It is well known and follows by a simple spectral argument, as Ah is self-adjoint 
positive semidefinite on S^, that the estimate 

(5.2) \\El k P h v - E h (t n )P h v\\ < CfciV'IMI, 7 6 [0, 1], 

holds for the backward Euler method P^ j . It follows from the stability of the finite 
element approximation and [11, Corollary 5.3] that 

(5.3) \\E h (t)P h v - E(t)v\\ < Ch-'t-^ 4 \\v\\, 7G[0,r]. 

Thus, with 7 = a/4 < r/4 < 1 in (jOj) and 7 = a in ([O]) . the estimate ([5~Tj) 
follows. 

It is also well known (see, for example, [HI Theorem 6.13]) that 

(5.4) \\(E(t)-E( S ))A-rv\\<\t- S \i\\v\\, 7 e[0,l], 
and therefore, taking also (|5.1[) into account, it follows that 

(5.5) \\(EhAt)Ph ~ E(t))v\\ < C(h a + k a / 4 )t- a / 4 \\vl a£[0,r\. 
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Indeed, for t G (tj-i,tj] we have that 

\\{E{t)P h - E{t))v\\ = \\(E{ k P h - E(t))v\\ 

< \\(E{ k P h - E{t 3 ))v\\ + ||(£fe) - E(t))v\\. 

For the first term (|5.1j) applies and for the second term we use (|5.4|) : 

\\(E( tj ) - E(t))v\\ = \\A a '*E(t)(E(ti - *) " ^""Nl 

< ||A a / 4 B(t)||||(B(t i -t)-J)i4- a / 4 t;|| < Ck a /H- a ' 4 \\v\\. 

Finally, we recall the smoothing property of the backward Euler scheme. It follows 
from |24[ Lemma 7.3] by stability and interpolation that for t G (tj—i,tj], 

\\A a h E Kk {t)P h v\\ = \\A a h Ei k P h v\\ < CtJ a \\Pkv\\ < Ct- a \\v\\, a>0. 
Therefore, 

(5.6) \\AlE Kk {t)P h v\\ < Ct- a \\v\\, a>0,t>0. 

We are now in the position to prove the following estimate for the weak error in 
case of the linearized CHC equation. As it was the case for the wave equation, the 
weak convergence rate is twice that of the strong convergence rate |17] (up to a 
logarithmic factor) under essentially the same regularity requirements on A and Q. 

Theorem 5.1. Let X be the solution of (|1.1[) and X£ k be given by (|4.23j) with 
spaces and operators described above and B^ = P^. Assume G G C^(H,M), Xq G 
Li(n,H) and 

(5.7) ||^- 2 >/ 2 Q|| Tr < K, \\A{ - 2)/2 P h Q\\ Tl < K, 

for some (3 G (0, £] and K > 0. Then there is C depending on T, K , ||-X"o||Li(fi,w)> 
and ||G|| C 2/ W j{-) such that, for Nk = T, h 4 + k < T, 

(5.8) |E(G(X&) - G(X(T)))\ < C(h^ + fc^ 2 ) log(^). 

Proof. Assumption (|5.7j) guarantees that || ^4 3 <5 ^ 1 1 hs < oo in view of Theorem 



12.11 This in its turn implies that X exists, as shown in [17 . Will use Theorem 
ETflwith O as in (|37f|) . Furthermore, let B = P h ,X = P h X , E(t) = E h>k (t) and 
Y htk (t) be defined as in (|4T24j) . whence Y h , k (T) = X* k . The use of Theorem (JS3J is 

justified since Eh,k(t)PhQ[Eh,k(t)Ph}* is bounded and of finite rank so that (|l.ll|l 
holds. 

We write F htk (t) = E h . k {t)P h - E(t) and recall (|ST5) . For the first term in ([376]) 
we use that Y htk (0) - F(0) = F htk (T), (|5"3|) . and the bound for u x in (|3"H) to get 

E j (u x (Y(0) + s(Y hik (0) - Y{0)),0),Y hik (0) - F(0)) ds 

( 5 - 9 ) < sup 11^(1, 0)|| E(||F M (T)X ||) 



xen 



< 



\ cUn , R) C(h^ + k^)T~^E(\\X \\ H ) 
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For the second term of (|3.6[) we have by (|2.3|) and repeated use of (|2.4[) that 
E J Tr (u xx (Y h , k (t),t)(E h>k (T - t)P h + E(T - t))QF Kk {T - t)*) At 

< E( J \\u xx (Y hik (t),t)(E htk (T - t)P h + E(T - t))Q\\ Tr \\F h ,k( T - *)lls(70 dt 

< sup \\u xx (x,t)\\ B ( H ) 

(i,t)E«x[0,T] 

x /" T ||(A-^- 2 >/ 2 ^, fc (i)^- 2 ^ 2 ^ + A-^- 2 )/ 2 ^(t)A^- 2 )/ 2 )Q|| Tr 
Jo 

x\\Fh,k(t)\\B(H)dt 

< \\G\\c*(u,m ^ (\\A^- 2)/2 E htk (t)Um \\A ( t 2)/2 P h Q\\Tr 
+ \\A-^- 2 ^ 2 E(t)\\ B{n) ||^- 2 )/ 2 Q|| Tr ) ||F M (t)|| B(70 dt 

(«, R )(ll^r 2)/2 ^QIlT r + ll^- 2 )/ 2 g|| Tr ) 

(||A-^- 2)/ Xfc(*)ll B («) + \\A-^-^ 2 E(t)\\ B{n) ) \\F htk (t)\\ m) At. 



< \\G\\a 



Jo 



\c?(h;. 



By (J5.7I) the factors in front of the integral are bounded by 2K\ 

We proceed by splitting the integral in two as L = L + J h 4 +k - For the first 
integral we notice that the last factor of the integrand is uniformly bounded and 
hence, by the analyticity of E(t) and (|5.6p with a — — ((3 — 2)/2, 



h*+k 




(\\A- h {ti - 2)/2 E Kk {t)\\ B(H) + \\A-^-^E{t)\\ B{n) ) \\F Kk {t)\\ B{n) At 



<C [ h+k t^-V' 2 At = C(h 4 + kfl 2 < C(h 2 ? + k^ 2 ). 
Jo 

For the second part we use again the analyticity of E(t), (|5.5[) with a — 2/3 and 
(|Ed| to get 



T 

h 4 + k 



(||^^- a )/ 2 £ M (t)|| B(w) + [[A-W- 2 )/ 2 £(*)|| B( «))||FM(*)llB(w)dt 
< C f {t^-V' 2 + i(^ 2 )/ 2 ) (h*f> + k^ 2 )t-^ 2 At 

Jh A +k 

;./3/2^ / 4.-1 A* — ri-\„„( 

h A + k ■ 

D 



C(h 2p + k^ 2 ) f r 1 At = C\og(j^)(h 2 P + k?/ 2 ) 

Jh^+k 



Remark 5.2. We refer to |15[ Theorem 4.4] for /i-independent conditions guaran- 
teeing (|5.7j) and the remarks after its proof for further discussions of the abstract 
conditions. Furthermore, the dependence on T of C in (|5.8|) can be removed if we 
assume that X$ £ Li(il, H 2 ^) by using the deterministic error estimate for smooth 
initial data from [T7| in (15.91). 
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Remark 5.3. The weak convergence of the finite element space discretization and 
backward Euler time discretization of stochastic heat equation with additive noise 
was considered in [TTJ]. The results there can be recovered using the fully discrete 
deterministic estimates 

\\E% k P h - E(t n )\\ B(n) < C(h 2 + k)t~ l 
and 

\\A a El k P h \\ B(n) + \\A a E(jtn)\\ B (H) < Ct- a , a e [0, §], 
together with Theorem 13.11 The technicalities are the same as in the spatially 
semidiscrete case [151 Theorem 4.1] under the same symmetric condition 

IIA^Q^Ihs < oo. 

We do not detail this here any further as it recovers a known result, only with 
perhaps a more transparent proof. 
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